A dynamic nonlinear model for saturation in industrial growth 
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A general nonlinear logistic equation has been proposed to model long-time saturation in industrial growth. 
An integral solution of this equation has been derived for any arbitrary degree of nonlinearity. A time scale for 
the onset of nonlinear saturation in industrial growth can be estimated from an equipartition condition between 
nonlinearity and purely exponential growth. Precise predictions can be made about the limiting values of the 
annual revenue and the human resource content that an industrial organisation may attain. These variables have 
also been modelled to set up an autonomous first-order dynamical system, whose equilibrium condition forms a 
stable node (an attractor state) in a related phase portrait. The theoretical model has received close support from 
all relevant data pertaining to the well-known global company, IBM. 
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In view of the current economic recession that is prevail- 
ing globally, it has become imperative to furnish mathemati- 
cal models of greater degrees of quantitative accuracy to un- 
derstand stagnation in economic growth, be it of states or of 
industrial organisations. When it comes to making forecasts 
about how the future might unfold for a state economy or an 
organisation, one realises the acute necessity for proper deter- 
ministic models. A recognition is gradually gaining currency 
that economic models should have robust predictive power, 
and this power should issue from the compatibility of the per- 
tinent mathematical models with empirical data 

Addressing this requirement is the principal objective of 
this work, and this has been done by making a case study 
on industrial growth. Many aspects of industrial growth lend 
themselves to well-formulated mathematical analyses. The 
health of a company is to be judged from the revenue that it is 
capable of generating, as well as the extent of human resource 
that it is capable of employing in achieving its objectives. One 
could make quantitative measures of all these variables, and 
this makes it relatively easy to have a clear understanding of 
industrial growth pattern, as well as to posit a mathematical 
model for it. The approach to these issues here is predom- 
inantly based on the use of standard mathematical tools of 
nonlinear dynamics and dynamical systems |12t|3t|^. There is 
a general appreciation that even when an industrial organisa- 
tion displays noticeable (very commonly exponential) growth 
in the early stages, there is a saturation of this growth towards 
a terminal end after the elapse of a certain scale of time |5|. 
As the system size (reflective of the scale of operations) be- 
gins to grow through the passage of time, a self-regulatory 
mechanism takes effective control over growth and gradually 
drives the system towards a saturated terminal state. And so, 
on very large scales of industrial operations, a clear under- 
standing could be derived about the constrained feature of the 
space within which an organisation functions. 

The theory developed here has been subjected to empirical 
test with the help of data collected from an industrial organi- 
sation that is global in character, i.e. its presence is to be seen 
and felt all over the world. This choice is dictated by the re- 
quirement that one would like to understand the global growth 



behaviour of a company, whose operating space is by defini- 
tion on the largest available scale, and, therefore, the overall 
pattern of its growth would be free of local inhomogeneities. 
To this end it has been worthwhile to study the revenue gener- 
ating capacity and the growth of the human resouce strength of 
the multi-national company, IBM. This organisation has been 
in existence in its presently known form for nearly a century. 
Besides this, it has spread all over the globe. So on both 
of these counts, a company like IBM is perfectly suited for 
the present study. Data about its annual revenue generation, 
the net annual earnings and the cumulative human resource 
strength, dating from the year 1914, have been published on 
the company website' itself . It has been satisfying to note 
that, analysed according to the stated objectives and speci- 
fications of this work, the IBM data actually give a striking 
match with the mathematical models forwarded here. Both 
the capacity for revenue generation and the human resource 
content of IBM, over a period of more than ninety years of the 
existence of the company (this long period is actually quite 
expedient for this study, since it is concerned with the growth 
of an industrial organisation from its inception to its termi- 
nal stage), show an initial phase of exponential growth, to be 
followed later by saturated growth towards a terminal state. 

A growth trend of this type can be described satisfactorily 
by a logistic differential equation, usually of second-degree 
nonlinearity, as it is done to study the growth of a popula- 
tion 10, lit]. Regarding the study of growth from industrial 
data, a preceding work |6] has pedagogically underlined the 
relevance of various model differential equations of increas- 
ing complexity. Along these lines, a generalisation of the lo- 
gistic prescription, to any arbitrary degree of nonlinearity, is 
being posited here, to follow industrial growth through time, 
t. Such a generalised logistic equation will read as 

</>(t) = A0(l-?7r) , (1) 

where can be any relevant variable to guage the health of a 
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firm (with the "dot" on (f> being its simple time derivative), Hke 
its annual revenue (or cumulative revenue growth) and human 
resource strength. The parameters a and 77 are, respectively, 
the nonlinear saturation exponent, and the "tuning" parame- 
ter for nonlinearity. Both influence the saturation behaviour 
of firm growth. In the context of the growth of an industrial 
organisation, one can identify various factors that may con- 
tribute collectively to this saturation in growth. These factors 
can be both economic and non-economic in nature. Some may 
operate internally, while others can make their impact exter- 
nally. A primary factor, regarding this study at least, is the 
space within which an organisation can be allowed to grow. If 
this space is constrained to be of a finite size (as, in a practical 
sense, it has to be), then, of course, terminal behaviour be- 
comes a distinct possibility. Even as it continues to grow, an 
organisation will gradually have to contend with the bound- 
aries of the space within which it has to operate. This brings 
growth to a slow halt. Indeed, saturation in growth due to 
finite-size effects is understood well by now in other situa- 
tions of economic interest where physical models can be ap- 
plied [7]. The adverse conditions against growth can be fur- 
ther aggravated by the presence of rival organisations com- 
peting for the same space. The last factor can become particu- 
larly crucial when a miscalculation is made in assessing future 
directions of growth vis-a-vis those of rival organisations — 
both the existing ones and the ones that might emerge in the 
future. 

Integration of Eq. ([T]i, which is a nonlinear differential 
equation, yields the general integral solution (for a ^ 0), 

4){t) = [?7 + exp i-aXt)] , (2) 

in which c is an integration constant. The fit of the forego- 
ing integral solution with the data has been shown in Fig. [T] 
which gives a log-log plot of the annual revenue, i?a, that IBM 
has generated over time, t. Here the annual revenue has been 
measured in millions of dollars, and time has been scaled in 
years. The data and the theoretical model given by Eq. (|2]i 
agree well with each other, especially on mature time scales, 
where nonlinear saturation is most conspicuous. Similar fea- 
tures are evident in Fig.|2] which plots the cumulative growth 
of the /BM revenue, i?c (in millions of dollars), with each suc- 
cessive year, plotted as t. The prescription of the cumulative 
growth needs to be stressed upon here. For any evolving sys- 
tem, its rate of growth is almost always a direct mathematical 
function of its current state (which has been built up in a cu- 
mulative fashion). Very frequently this kind of dependence 
leads to an exponential growth pattern, something that is quite 
relevant here, especially as regards the early growth of IBM. 
Besides, the distribution of the cumulative revenue is much 
more free of fluctuations than the distribution of the annual 
revenue, and as such the former is better suited for modelling. 

While the exponential feature may be appropriate for mod- 
elling the early stages of growth, the later stages shift into a 
saturation mode. A limiting value in relation to this saturated 
state could be found by setting (/) = in Eq. ([T]), and this will 




FIG. 1: The continuous curve gives the model fit for the annual rev- 
enue, _Ra, generated by IBM. The fit by the theoretical model agrees 
well on nonlinear time scales, for a = 1, A = 0.145 and 77 = 10~^. 
Here 7?a has been scaled in millions of dollars, while t has been 
scaled in years. 




FIG. 2: The cumulative growth of the annual revenue, R^, generated 
by IBM is fitted by the theoretical model (given as the continuous 
curve). The fit is very close once again on nonlinear scales, for a — 
1, A = 0.15 and rj — 5 x 10~^. As before, Rc has been measured in 
millions of dollars, while t has been scaled in years. This distribution 
is evidently more free of fluctuations than the previous one. 



lead to 

0sat = 77"'/" • (3) 

Making use of the values of a and 77, which have fitted the 
saturation properties of the plot in Fig.[T] a prediction can be 
made that the maximum possible annual revenue {not the cu- 
mulative revenue plotted in Fig.|2]i that IBM can generate will 
be about 100 billion dollars. 

Another point of great interest is that the growth data have 
been fitted very well by the simplest possible case of non- 
linearity, given by a = 1. This will immediately place the 
present mathematical problem in the same class of the logistic 
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FIG. 3: The net annual earnings made by IBM has shown steady 
growth, except for the early years of the 1990s decade, which was 
about 80 years of the existence of the company. Around this time the 
company suffered major losses in its net earnings, and this time scale 
corresponds very closely to the time scale for the onset of nonlinear 
saturation in growth indicated by Eq. l|4). 



differential equation devised by Verhulst to study population 
dynamics [2,^]. This equation has also been applied satisfac- 
torily to a wide range of other cases involving growth • That 
the industrial growth model, based on the IBM data, should 
also fall within the same category is worthy of special note. 

The time scale for the onset of nonlinearity can be deter- 
mined by requiring the two terms on the right hand side of 
Eq. (|2|i to be in rough equipartition with each other. This will 
give the nonlinear time scale, 

t„i -(aA)-Mnhc"|, (4) 

from which, making use of the values of a, 77 and c needed 
to calibrate the IBM revenue data (both the annual revenue 
and the cumulative revenue), one gets t^i 75 — 80 years. 
It is easy to see that this is the same time scale on which 
nonlinearity causes the onset of saturation in the growth of 
the company, and an indirect confirmation about the valid- 
ity of this time scale comes from the plot in Fig. |3] which 
shows the growth of the net annual earnings of IBM (labelled 
as P, the profit), against time, t. The company suffered ma- 
jor reverses in its net earnings (upto 8 billion dollars in 1993) 
around 1991-1993, which was indeed very close to 80 years of 
the company, since its inception in 1914. This intriguing cor- 
respondence between the two time scales, arrived at via two 
distinctly different paths, is arguably much more than a simple 
coincidence. 

The human resource content of an industrial organisation 
is also another indicator of its prevailing state. In the case of 
IBM, the data for the human resource of the company have 
been plotted in Fig. H] Over the years the growth of human 
resource has been steady, and has followed the qualitative 
growth trend of the revenue. However, one may readily no- 
tice that just like the net earnings of IBM, there has been a 
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FIG. 4: The growth of the human resource strength is fitted globally 
by the theoretical model for a = 1, X — 0.09 and rj = 2 x 10~®. 
There has been a noticeable depletion of human resource on the same 
nonlinear saturation time scale given by Eq. i.e. 75 — 80 years. 



sharp decline in the number of human resource around the 
time scale of 75 — 80 years of the age of the company. This 
may be construed as another independent piece of evidence in 
support of the theoretical estimate of the saturation time scale 
that Eq. (HI indicates. And going by the values of a and ij 
needed for the model fit in Fig. 21 the maximum possible hu- 
man resource strength that IBM can viably employ is predicted 
from Eq. © to be about 500, 000. 

By now it has become evident that there is a strong cor- 
relation among the variables by which one may monitor the 
state of an industrial organisation. And in fact it happens that 
a complete impression of the true state of affairs may only be 
conveyed by analysing all the variables involved. The concept 
of the "Balanced Scorecard" is somewhat related to this prin- 
ciple (si]. The meaningful variables can all be diverse in na- 
ture, ranging from the finances of an organisation to its human 
and technological resources and to the market within which 
the organisation might be operating. The growth rate of any 
one of these variables may have a correlated functional de- 
pendence on the collective current state of all the pertinent 
variables. 

So if one were to study an industrial organisation whose 
current state is defined completely by a general set of n vari- 
ables, like Rc (or R^) and H, then the growth rate of the i-th. 
variable, 4>i, will be given by 

<^>^ = ^^iWJ}), (5) 

with being a general function of all the variables in the set. 
The whole set can be expressed explicitly by making both i 
and j run from 1 to n. This will give a set of n first-order 
differential equations, with each one of them being coupled 
to all the others, and this entire set will form an autonomous 
first-order dynamical system in an n-dimensional space. 

If an industrial organisation generates enough revenue, it 
becomes financially viable for it to maintain a sizeable human 
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resource pool. On the other hand, the human resource strength 
will translate into a greater ability to generate revenue. In this 
manner both the revenue and the human resource content of 
an organisation will sustain the growth of each other. Defining 
a general revenue variable, R (which can be either i?a or Rc), 
its coupled dynamic growth features along with the human 
resource, H, can be formally stated in mathematical terms as 



R = p{R,H) 
H = a{R,H) 



(6) 



The foregoing coupled set of autonomous first-order differen- 
tial equations forms a two-dimensional system involving two 
variables, R and H. The equilibrium condition of this dynam- 
ical system is obtained when both the derivatives on the left 
hand side of Eqs. (|6]l vanish simultaneously, i.e. R — H — 0. 
The corresponding coordinates in the H — R plane may be la- 
belled {Hq, Rq). Since the terminal state implies the cessation 
of all growth in time (i.e. all derivatives with respect to time 
will vanish), it is now possible to argue that the equilibrium 
state in the H — R plane actually represents a terminal state in 
real time growth. 

Some general deductions can now be made about the nature 
of the equilibrium state, with the help of dynamical systems 
theory p, '7). The two coupled equations, given by Eqs. 
will, in the most general sense, be nonlinear. A linearisa- 
tion treatment on them could be carried out by applying small 
perturbations on R and H about their equilibrium state val- 
ues. The perturbation scheme will be i? = Rq + R' and 
H = Hq + H' . This will allow a coupled set of linearised 
equations to be set down as 



AR' 
dH' 



AR' + BH' 
CR' + VH', 



(7) 



in which 
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(8) 



Solutions of the form R' 



and H' 



will enable one 



to derive the eigenvalues, uj, of the stability matrix implied by 
Eqs. dTji, as 



{A + V)± V {A + P) - 4 (AV - BC) 



(9) 



The exact determination of the values of the two roots of lo 
will impart a clear idea about the nature of the equilibrium 
state, which can either be a saddle point or a node or a fo- 
cus |l3l|4|]. The last case will necessarily mean an oscillatory 
nature in the growth of both R and H through time and 
is not tenable here. Going back to Figs.[ri|2]&|4] one notices 




FIG. 5: The straight-line fit validates the logistic equation model. 
The slope of the straight line is given as 1.4, and it closely matches 
the value of 1.6 that could be obtained from the parameter fitting in 
Figs.|2]&|4] The cusp at the bottom left is due to the loss of human 
resource. The lower arm of the cusp has nearly the same slope as the 
straight-line above it. Growth, corresponding to the positive slope in 
this plot, can be modelled well by the logistic equation. 



that the respective growth patterns of both R and H have, on 
the other hand, been largely monotonic in nature (except for a 
manifestly sharp dip in H at large values of t). So an immedi- 
ate conclusion that follows is that the equilibrium state is very 
likely a node flHl, and this will correspond mathematically 
to u) having two real roots of the same sign. Practically speak- 
ing, this is what is to be expected entirely. For an industrial 
organisation it is not conceivable that while there is growth in 
one variable, there will be decay in the other (if this were to 
happen, the equilibrium state will be a saddle point). Both will 
have to grow in close mutual association, and, if the dynam- 
ical systems argument is anything to go by, both will make 
the approach towards the terminal state simultaneously. The 
qualitative behaviour of one variable cannot be completely in- 
dependent of the other. 

To have a quantitative corroboration of the dynamical sys- 
tems modelling with the help of the IBM data, one could, mo- 
tivated by Eq. ([T]i, proceed with a simple ansatz to describe 
the individual growth patterns of R and H by an uncoupled 
logistic equation model. This can be set down as 



R{t) = XrRil-VrR"'-) 

Hit) =. XhH {1 - r^hH"^) 



(10) 



with the subscripts r and h in the parameters a, A and rj, indi- 
cating that R and H will each, in general, have its own differ- 
ent set of parameter values. From Eqs. ( fTOl i one could easily 
eliminate the two time derivatives and arrive at 



dR 
dH 



\rR{l-T]rR°'- 



(11) 



XhH {1 - r]hH^<^) 
whose integral solution can be expressed in a compact power- 
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law form as 



under the definitions that 



V = Vr ■ 



(12) 



, (13) 



and n is an integration constant. The power law in Eq. (fTZt 
implies that a log-log plot of v against u will be a straight 
line with a slope, (3. For the greater part of it, this fact has 
been established emphatically in Fig. |5] In this plot, v has 
been defined in terms of the cumulative revenue, i.e. R = Rc, 
while going by the model fitting in Figs.[T]|2]&|4l the values 
ar — ah — 1, rjr = 5 X 10^'' and rjh = 2 x 10^^ have 
been used. The cusp in the bottom left corner of the plot has 
arisen because of an irregular depletion of human resource in 
IBM in the early 1990s, and hence it is beyond the scope of 
the logistic equation model. However, the lower arm of the 
cusp has nearly the same positive slope as the straight-line fit. 
So, as long as there is steady growth in both R and H, the 
logistic equation model gives a good description of the global 
performance of an industrial organisation. A company can be 
said to be performing well if it is on a curve with a positive 
slope in the u — v plot, as it has been implied in Fig.|5] Sus- 
tained deviations from the straight line and having a negative 
slope in this log-log plot, are causes for worry regarding the 
well-being of the company. The slope of the straight line in 
Fig. |5] sets a value for the power-law exponent in Eq. (fT2l l as 
(3 ~ 1.4, which is quite close to the value of /3 ~ 1.6, found 
simply by taking the ratio of the respective theoretical values 
of A, chosen to fit the empirical data in Figs.|2]&|4] 

It is remarkable that the simple logistic equation employed 
for the dynamical system in Eqs. ( fTOb . would suffice to model 
the IBM data so closely, as Fig.|5]shows in no uncertain a man- 
ner Once it is obvious that the dynamical systems modelling 
has been successful, one could formulate a precise description 
of the terminal state. This can be done by referring to Eqs. (O 
and (|9]l, for which, going by Eqs. (fTOl i. it is easy to see that 
neither does p depend on H, nor does a depend on R. Con- 
sequently, one has B = C = Q. Immediately the two roots 
of Eq. (|9]l are found to be w = A and uj ~ T>. To ascertain 
the exact values of A and V in terms of the known parame- 
ters, one will have to know the equilibrium values of R and 
H. These can be shown to be Rq — -q^^ and = r/^^, 
corresponding Xo u — v — Q for the terminal state. After this, 
it is straightforward to argue that A — —\r and B = — A^. 
Since the values of A are positive in all the cases of model fit- 
ting, it is now possible to claim that with both the roots of u) 
being real negative numbers, the limiting state for industrial 
growth is represented by a stable node in the phase portrait of 
an autonomous first-order dynamical system Extending 
this contention further, the limiting state can be perceived to 
be an attractor state, towards which there will be an asymp- 
totic approach through an infinite passage of time ['T). At least 
this is what the growth data pertaining to IBM indicate. 



The terminal feature in industrial growth derives from the 
nonlinear term in the general logistic equation used for the 
modelling. This will naturally suggest that the retarding fac- 
tors acting against the growth of an organisation are nonlinear 
in character. These factors are all quantified through the pa- 
rameters a and -q in Eq. ([T]|. Of these two parameters, 77 is 
more amenable to manipulations, and conditions conducive 
to lasting growth will necessitate tuning it down. This ought 
to be the guiding principle behind a successful management 
strategy for feasible long-term growth, especially in the case 
of organisations that are still in their early stages. Knowledge 
of the general nature of the possible adversities lying ahead, 
can enable a company to apply appropriately corrective mea- 
sures at the right juncture, and defer the arrival of the nonlin- 
ear time scale. As a result it will make a more effective im- 
plementation of future strategies and innovative solutions for 
growth, all of which should be in a state of adaptable align- 
ment with their objectives and core competencies. The "Blue 
Ocean" strategy [9], for instance, might be one such practical 
and effective means. 

On the other hand, it would never be too late for companies 
which would already have entered a terminal phase, to insti- 
tute proper reality checks. This will allow for a more func- 
tional and timely redefining of fundamental objectives. The 
solutions which might follow, could be varied in many un- 
expected ways. Rather than adhering to conventional modes 
of growth and survival, industries could devise ways of pre- 
serving both their existence and their relevance by being more 
integrated with the social welfare of their markets — in short, 
assist in creating an environment of overall prosperity and a 
general feeling of well-being, right alongside the creation of 
wealth. The benefits derived through such inclusive strategy 
implementation should be lasting. 
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